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ABSTRACT 


This thesis contains a study of the interannual and intra- 
annual tropical atmospheric circulations for the Northern 
Hemisphere summer seasons of 1974, 1975 and 1976 using satel- 
lite observations of infrared radiation and albedo in conjunc- 
tion with parameters available from conventional 200 mb data. 
The 200 mb flow regimes are contrasted for the three seasons 
and power spectra and cross-spectral analysis are used to 
find significant correlations between the parameters from 
time series generated within a given area and between the 
areas. These areas were picked for each parameter to represent 
regions of maximum time-averaged values. The anomalously warm 
sea Pesace temperatures in 1976 (Bepristis, 1977) occurred 
in conjunction with major circulation changes paralleling 
Bjerknes' (1969) hypothesis and the differences in flow re- 
gimes between a normal and dry monsoon season are in general 
agreement with Kanamitsu and Krishnamurti (1978). Signifi- 
cant correlations were found to exist between parameters with- 
in the eastern Pacifie and western Pacific and also between 
these two regions. Particularly strong correlations involved 


the i, vi and IR parameters. 
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i LMRODUCT ION 


Previous studies (Chang and Miller, 1977; Kanamitsu 
and Krishnamurti, 1978; Winston and Krueger, 1977; Bepristis, 
1977) have analyzed upper air and satellite data in an 
effort to study the interannual variations of planetary- 
scale circulation over the tropics. 

Chang and Miller (1977), in a spectral analysis of 
tropical easterly waves during the second half of 1972 and 
iews, found considerable variation in the structure and 
properties of the disturbances between the two seasons. In 
the tropical eastern Pacific the 1972 season was character- 
ized by abnormally large warm sea temperature anomalies, 
which had an effect of almost completely removing the nor- 
mal zonal sea surface temperature gradient, while the 1973 
season showed abnormally large cold sea surface temperature 
anomalies, which sharpened the gradient in the central 
Pacific. The periodicity of the waves in both seasons was 
found to be 4-5 days. However, the vertical structure and 
amplitudes of the waves exhibited marked differences. 

An extensive study completed by Kanamitsu and Krishna- 
fee (1978) contrasted the 200 mb flow regimes during a 
drought year over central India and western Africa (1972) 
with those during a normal rainfall year over the same region 


(1967). Notable differences during the drought year (1972) 


ING 





included the warm sea surface temperatures (Chang and Miller, 
1977) which they related to a near doubling of typhoon days 
over the western Pacific and a shifting of the Tibetan High 
by some 10° eastward and equatorward. 

Another study by Bepristis (1977) examined radiosonde 
@emea tor 1974, 1975 and 1976 for the islands of Majuro, 
Ponape, Truk, Yap, Koror and in the tropical central and 
western Pacific. He also constructed sea surface tempera- 
mme@e (Sol) anomaly charts for the tropical Pacific (Figure 
19) showing 1974 to be near normal, 1975 to be cold and 1976 
to be warm. He found the thermal structure of the easterly 
waves was influenced more by the immediate upstream (about 
mo longitude) SST than by the local SST. 

A fourth study, that of Winston and Krueger (1977), in- 
troduced global patterns of outgoing long-wave radiation, 
albedo, absorbed solar radiation and net radiation derived 
from satellite observations. These patterns along with 
patterns of interannual differences in seasonal heating were 
examined for the three summers 1974-76 over the Eastern 
Hemisphere in relation to the variations in the summer mon- 
soon and the major circulation features. Relating low 
values of outgoing long-wave gad ienen to cloudiness and 
rain they concluded 1975 was the strongest monsoon season, 
Sepecially over the south and east. In 1974 the lowest 
outgoing radiation occurred in the area of NE India, Nepal 
and the Himalayas resulting in decreased cloudiness and 


Meee all Over southern India. The 1976 monsoon had a 


eal 





delayed start, was fairly active during July and August and 
receded somewhat early in September. A possible link rela- 
tive to the delay of the monsoon in 1976 (and consequently 
weaker) compared to 1975 showed up in the equatorial Pacific 
eastward of New Guinea. In 1976 this area was characterized 
by decreases in long-wave radiation indicating more cloudi- 
ness and convection. Thus, the weaker monsoon was accom- 
panied by major changes in the regions of condensation heat- 
ing over the tropics. As we have already seen, (Kanamitsu 
and Krishnamurti, 1978), the increased cloudiness and convec- 
tion in the eastern Pacific should indicate warmer sea 
surface temperatures in 1976 and this was verified (Bepristis, 
may 7). 
The main objective of this work is to combine some of 

‘ 
the parameters available from satellite data (Winston and 
Krueger, 1977) with other parameters available from 200 mb 
data (Kanamitsu and Krishnamurti, 1978) for the summers of 
1974-76 to study interannual and intraannual fluctuations in 
the parameters ee qe divergence, relative vorticity, infra- 
red radiation and albedo. These parameters will be compared 
statistically to determine their inter-relationships within 


a particular system as well as between geographically sepa- 


rated systems. 
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Li eee A Dee NAbY SiS TECHNIQUE 


A. DATA SOURCES 
1. Satellite Data 
Satellite data were obtained from the National En- 
vironmental Satellite Service/NOAA, which provided digitized 
information of infrared radiation (IR), available solar 
radiation Soa? and absorbed solar radiation (S34? for the 
entire globe. The data are in radiative flux form, with 
values in watts per square meter. Area-averaged values were 
available on a 125x125 hemispheric grid with each grid value 
representing an area of approximately 150 km x 150 km. 
2. Upper Air Data 
200 mb U-component and V-component data from analyzed 
wind fields were obtained from the Fleet Numerical Weather 
Center (FNWC), on their 49x144 mercator global band grid. 
This grid covers the area from ud° south to 60° north. The 


grid mesh distance at 22.5° north/south is 257 km. 


B. DATA PERIOD 

All data were collected at 12 hour intervals for the 
Memtod 16 May through 15 September 1974, 1975 and 1976 with 
the exception of Sap and Se which were only available once 
per day and IR, Sab and =e for the period 15-31 May 1974 


Meer nOt available). 


digs: 





C. SRARAMEDRERS CAGCULATED 
Seasonal mean fields at 200 mb for the following para- 
meters were produced for the period 16 May - 15 September 
for each of the years 1974, 1975 and 1976. 
1. U-component 
Seasonal mean U-component fields were produced on 
a 49x144 mercator global band aes Simply summing up 
values at each grid point for each time period and dividing 
by the number of time periods. See Figures 1-3 for 1974, 
1975 and 1976 respectively. 
2. V-component 
Seasonal mean V-component fields were produced on 
a 49x1l44 mercator global band grid by simply summing up 
values at each grid point for each time period and dividing 
by the number of time periods. See Figures 4-6 for 1974, 
1975 and 1976 respectively. 
3. Divergence 
Seasonal mean divergence fields were calculated on 
a 49x1l44 mercator global band grid from the seasonal mean 
U-component and seasonal mean V-component fields using equa- 
ion (1) (Krishnamurti, 1970) with centered finite differenc- 


ing. See Figures 7-9 for 1974, 1975 and 1976 respectively. 


; oil av vy tan ¢ 
Divergence = > + 7 (1) 


Mimjree 2 1s the radius of the earth (=6371.229 km), (') is 


the time-averaged mean and @ the latitude. 
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4, Relative Vorticity 
Seasonal mean relative vorticity fields were calcu- 
lated on a 49x144 mercator global band grid from the 
seasonal mean U-component and seasonal mean V-component 
fields using equation (2) (Krishnamurti, 1970) with centered 


memete differencing. 


du _ Gl seuel w) (2) 


Relative vorticity = EF 5 


ov 
ax 
See Figures 10-12 for 1974, 1975 and 1976 respectively. 
ee ocream Function 

Seasonal mean stream function fields were calculated 
from the seasonal mean relative vorticity fields by solving 
the Poisson equation (3) (Krishnamurti, 1971), with boundary 
conditions given by equation (4) (Brown and Neilon, 1961). 


See Figures 13-15 for 1974, 1975 and 1976, respectively. 


V°p = seasonal mean relative vorticity G2) 
Ve 2 ie (4) 
ae oF j=1,2,---,N 
where ae velocity component normal to the boundary 
A = starting value minus value at the end of integra- 


tion around the perimeter 


[Romnumoc Of CiStinect ead points around the 
boundary 


Ins 





6. Velocity Potential 
Seasonal mean velocity potential fields were calcu- 
lated from the seasonal mean divergence fields by solving 


mem romsson equation (5)"(Krishnamurti, 1971). 


VX = seasonal mean divergence (5) 


xy was assumed = 0 on the northern and southern boundaries 
Saemene grid and cyclic boundary conditions were used along 
the east-west side. See Figures 16-18 for 1974, 1975 and 


1976 respectively. 


Meee AREAS OF INTEREST 

Bao ecocnw ot cme 200 mb sseasondl mean U-component, V- 
component, divergence and relative vorticity fields, (Figures 
ieee) three areas of interest were selected. From the 
divergence fields, the three areas were picked to represent 
regions of maximum time-averaged divergence and from the 
relative vorticity fields the three areas included the maxi- 
mum negative vorticity centers. The three areas of interest 
on the U-component and V-component fields were picked as 
Mme recions of maximum North-South gradient of isopleths of 
stream function and velocity potential respectively. As a 
result of the selection criteria, areas may be different for 
different parameters and different seasons. Also,different 
areas may relate to the same system, e.g. easterly jets, 


anticyclones, western or eastern Pacific ITCZ, etc. 
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Peete SERIES PRODUCED 

With the areas of interest defined and a time period of 
16 May - 15 September in 12 hour intervals specified, 18 
time series of area-averaged values were produced for each 
year. Three each from the U-component and V-component data 
where time series of the parameters we and tie respectively 
were calculated instead of U and V as being more related to 
kinetic energy of the flow and their magnitudes less affected 
by the area-averaging process. There also were three of 
divergence and three of relative vorticity. The final six 
were from the digitized satellite data and were calculated 
over the same areas used for the divergence. Three of these 
were of area-averaged IR and three of area-averaged albedo 


using 


albedo = —24 ab 


av 
where Say 1s the available solar energy and 5 1s the ab- 
sorbed solar energy. A nine-point low pass filter was 
applied to all time series to remove the seasonal trends 


and eliminate periods in excess of 20 days. 


meee SLATLSTICAL ROUTINE 

The filtered time series were input to one of the routines 
of UCLA's Biomedical Computer Programs (BMD) Package (Dixon, 
mee, BMDOZT. Each time series was cross-correlated with 
all others. Parameters output included power spectrum esti- 
mates, coO-spectra, coherence squares and cross covariances. 
A lag of 40 was chosen providing resolution of periods up to 


20 days. 
Hy 








ie weteel lal S 


Peco rARISGON OF THE THREE YEARS 

In this section we shall examine the major differences 
in the large-scale 200 mb circulations during the years 
1974-76. Winston and Krueger (1977) have already classified 
the three years relative to the Summer Asian Monsoon region 


as being dry (1974), normal (1975) and late starting/weaker 


mee S). 
ime, Mean Circulactons 
Figures 1-3 show the fields of u for the years 1974- 
76. The intensity of the tropical easterly jet over the 


Arabian Sea is nearly the same for all three years (around 
eo ml a) but the jet core is displaced about 6° north and 
7° west in 1974 from the nearly coincident positions of 1975 
and 1976. Winston and Krueger (1977) observed a similar 
displacement of the maximum cloudiness and precipitation 
areas in 1974. Kanamitsu and Krishnamurti (1978) observed 
a Significant decrease in the intensity of the tropical 
easterly jet during a drought year over the monsoon region 
mes 72). 

Mid-latitude westerlies around 40° north are weakest 
in 1974 becoming stronger in 1975 and strongest as well as 
Shifted eastward in 1976. This is consistent with an inten- 


Sification and a shifting eastward of the Hadley-type 


Is 





| 
| 
| 





Circulation in 1976 as indicated by the locations of the 
positive centers in the velocity potential (Figures 13-15), 
stronger 200 mb divergence (Figures 7-9) and increased pole- 
ward motion over the eastern Pacific. Also consistent with 
an eastward shift in the jet core is the weakening of the 
memor 200 mb anticyellone in 1976 as indicated by the mean 
stream function (Figures 16-18). 

In 1975, westerlies around 20° north over the mid 
Pacific are weakest and the equatorial easterlies are broken 
in the region 140°-175° west longitude. This 1s consistent 
with the southwest penetration of the mid Pacific trough 
(Figure 17) indicating upper-level cyclonic circulation in 
the mid Pacific. This shows up on the mean divergnece field 
(Figure 8) as a strong divergent center and results in an 
intensified equatorial Walker circulation as seen from the 
mean velocity potential field (Figure 14). By contrast, 
1976 has the strongest westerlies around 20° north in the 
mid Pacific. Looking at the stream function (Figure 18) we 
See ridging across the entire Pacific moving the mid Pacific 
feeeueh to a poSition north and east of 1975. The appearance 
See SLYONg poSitive vorticity center (Figure 12) over this 
mid Pacific region indicates increased anticyclonic circu- 
lation, decreased divergence (Figure 9) and a strengthening 
Seeene Casterlies forming a continuous band across the equa- 
(jerial Pacific. This supports Kanamitsu (Kanamitsu and 
Krishnamurti, 1978) who related strong westerlies in this 
region to a drought year when comparing the years 1967 and 


ie 2. 


i 








The region of poleward motion in the eastern Paci- 
fic (Figures 4-6) extends further west and is intensified 
in 1976 consistent with warmer sea surface temperatures and 
increased convection, latent heat release implied from the 
stronger upper-level divergence (Figures 7-9) found over 
these regions. In the eastern Mediterranean region, the 
poleward motion center is shifted southeastward and inten- 
sified in agreement with the southeastward movement of the 
main 200 mb anticyclone, the decrease in cloudiness and 
precipitation and a resulting increase in solar insolation. 

2. Divergence 

Figures 7-9 illustrate divergence for the years 1974- 
76. The strongest year 1976 has a nearly continuous band of 
divergence across the Pacific between 5° and 10° north in- 
‘dicating a much more active ITCZ probably as a result of 
the increased sea surface temperatures (Figure 19). This 
Supports Winston and Krueger's (1977) areas of reduced long- 
wave radiation in 1976. The strong divergent center in the 
mid Pacific for 1975 reflects inereased convection, latent 
heat release and stronger cross equatorial circulation as 
seen in the velocity potential (Figure 14). The divergence 
center associated with the major anticyclone is shifted 
northwest and is more intense in 1974. This agrees with the 
observed shift in cloud and precipitation patterns (Winston 
and Krueger, 1977) and reflects the increased insolation 


over a dry southern India. 
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24) Relaweve Vomticity 

Figures 10-12 illustrate relative vorticity fields 
for the years 1974-76. The mid Pacific in 1976 shows a 
stronger area of relative vorticity indicating an increased 
anticyclonic circulation in the region of the observed in- 
crease in the easterlies. The strongest negative relative 
vorticity in all three years falls just north of a line 
connecting the vortices of the major anticyclone in the 
stream function (Figures 16-18). 

mee velocity Potential 

Figures 13-15 show velocity potential for the years 
1974-76. We can see the increased Hadley-type circulation 
in 1976 in the western Pacific resulting in an increase in 
the mid latitude westerlies at the northern extent of the 
Hadley cell. Note the shifting of this main divergent center 
eastward and southward in 1976 with an eiongated center. 
Such a change perhaps reflects to a certain degree the warm 
sea surface temperatures across the Pacific and the increase 
in convective heat sources indicated by increased divergence 
to the east. We also see the east-west overturning between 
the main divergent center and eastern Africa is strongest in 
1976. This indicates increased subsidence over eastern 
Africa resulting in decreased cloudiness and precipitation 
as observed by Winston and Krueger (1977). The changes in 
positions of the positive velocity potential centers can be 
Mieectpreted as a shift eastward and intensification of the 


rising branch of the equatorial Walker cell. We also make 


Za 





note of weaker westerlies in the mid latitudes in 1975. The 
convergent center at eos a conedatent feature of 
all three years indicating cross equatorial flow. 
Se ond tvemgenin CUreculatlons 

The rotational part of the wind is computed from the 
stream function field. Figures 16-18 show the stream func- 
tion field for the years 1974-76. The anticyclone over the 
monsoon region has a split center in all three years. This 
anticyclone is weakest in 1976 consistent with the shift 
eastward of the maximum westerlies to correspond to the 
northern extent of the increased Hadley circulation and de- 
creased cloudiness over the Indian subcontinent. Along with 
the stronger mid latitude westerlies in 1976 there is ridging 
from the anticyclone nearly all the way across the Pacific. 
Mis ridging forces the mid Pacific trough to a position 
north and east of its 1975 location increasing the anticy- 
Glonic circulation in the mid Pacific and allowing the 


easterlies to dominate the equatorial Pacific. 


Pee FOWER SPECTRA 

The power spectral estimates are tabulated in Table I 
which shows the number of areas of each parameter that had 
a peak period within the range shown at the left. Thus, 
the number 3 indicates all three areas selected for that 
parameter in that year showed a peak period within the range 
Meeeed and a 2 indicates that 2 of the 3 areas selected had 


meeemcdency for a peak period within the range listed, etc. 


cae 
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indicates predominant period was 
4-5 days 


Table I. Power Spectral Estimates 
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There was a strong diurnal cycle in the IR but periods less 
than 3 days have not been included. Note the predominance 

of the 3-4 day periodicity in 1974 and 1976. The predominant 
period in 1975, considered the normal monsoon season, was in 
the 4-5 day range. A secondary period in the 5-7 day range 
Bisa snows Up well in 1974 and 1976 while the trend in 1975 
appears to have been toward a mean 4-5 day period. In 1976, 
the most active eastern Pacific year, there were more 
occurrences of a longer (greater than 8 day) period. MThis 
season was characterized by a more active ITCZ region and 
warmer sea surface temperatures across the Pacific equatorial 
regions (Bepristis, 1977). Sample plots of power spectral 
estimates are shown in Figures 20-23. Figures 20 and 21 

are for the maximum time-averaged divergence areas in the 
eastern and western Pacific (1976) respectively. Figures 22 
and 23 are for the maximum time-averaged ic in the eastern 
and western Pacific (1976) respectively. Figures 20 and 21 
show a maximum peak in the 3 to 4 day range and Figure 20 has 
a secondary peak in the 5 to 7 day period. Figure 22 has a 
maximum peak in the 5 to 7 day range with secondary peaks 

1 the 3-4 day and greater than 8 day range while Figure 23 


has a 3-4 day period. 


SCOURS RENCE SQUARES 
In this section we will be concerned with three things. 
First, we will summarize the relationships of the parameters 


whose areas fall within the monsoon system. Secondly, we 
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will look at the relationship of areas within the mid to 
eastern Pacific and finally the relationships existing be- 
tween the mid to eastern Pacific and the monsoon areas. 
Results are tabulated in Tables II-IV for the monsoon sys- 
tem (74-76), the mid to eastern Pacific system (74-76) and 
relationships between the monsoon and mid to eastern Paci- 
fic systems respectively. An entry in the table indicates 
the two parameters are correlated with at least a 95% 
level of confidence. The period at which this correlation 
exists is given on the left of the table in days and the 
phase of the transfer function given in parts of a circle 
mmpetne entry 1tself. From the table it is obvious we do 
fee GOonsider correlations for periods less than 4 days (or 
Meme than 20 by virtue of the filter used). 
1. Monsoon System Relationships 

Within the monsoon region, 1974 shows the best cor- 
relation of the parameters, mostly at the 99% confidence 
level. The divergence corresponded well with the u- para- 
meter as well as with IR and albedo. The ve parameter showed 
a strong correlation to i amnaualbedeo.) Site dry year 1976 
Was next in terms of correlation of parameters, mostly at 
the 95% confidence level. Here, the divergence showed a 
Beed correlation to the ie IR and albedo parameters. The 
MOrmal year, 19753, showed a correlation of divergence to 


2 tO, albedo. The area of 


IR and albedo parameters and in U 
Betative vorticity that fell within this region in each of 


the three years showed no correlation to other parameters 
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Pee Nor greater, level. Other than the close relation- 
ship of infrared radiation, albedo and divergence which 
were all calculated over the same areas, there does not 
appear to be any pattern to the correlations between the 
parameters or the periods within the parameter. 
2. Mid to Eastern Pacific System Relationships 

In this region, 1976 showed the best results in cor- 
relation of the parameters. This was probably to be expected 
since this was the strongest and most active year in the 
eastern Pacific. Once again the divergence showed a strong 
relationship to the IR and albedo parameters and also to na 
In 1974, only two parameters had areas of interest within 
this region (relative vorticity and v*) and they showed a 
meeawcorrelation. Surprisingly, in 1975 there were five 
parameters with areas of interest within this region but 
there were no correlations among them at the 95% level. 
Comparing the periods and phase relationships of the para- 
meters in Table III that are also in Table II we see the 
periods are the same but the mid to eastern Pacific system 
is generally more out of phase. Year to year comparisons 
of periods and phase in the mid to eastern Pacific are not 
possible due to the lack of correlations in the years 1974 
ome 1975 . 

3. Inter-System Relationships 

Because nearly all areas of interest chosen in 1974 

fall within the monsoon region, it is impossible to obtain 


Mier —system relationships for that year. Significant 
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correlations between the mid to eastern Pacific region and 
the monsoon system exist in both 1975 and 1976. In each 
year the v2 parameter in the mid to eastern Pacific corre- 
lates well with the IR parameter in the monsoon region but 
the period and the phase are different. The normal season 
1975 has the relative vorticity in the mid Pacific correlat- 
ing well with vom IR and relative vorticity in the vicinity 
of the monsoon area. During the dry season 1976, we have 
the relative vorticity in the mid Pacific related to the be 
and albedo parameters in the monsoon region and IR in the 
mid Pacific related to the U parameter of the monsoon region. 
While we were able to show some correlations, they 
certainly are not conclusive since the area selection 
criteria we used did not guarantee that we would be calcu- 
lating all parameters in both the mid Pacific and the 
monsoon regions. We feel more meaningful results can be 
obtained in future studies with more attention given to 


area selection. 
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TV. CONCLUSIONS 


In this study of the interannual variations of the 
tropical atmospheric circulations we have seen that in an 
anomalously warm sea surface temperature year a stronger 
Hadley type circulation developed and was shifted eastward. 
This was consistent with increased convection associated 
with a more active ITCZ, increased latent heat release, 
increased divergence aloft, a northward transport of 
momentum with a corresponding increase in the mid latitude 
westerlies and increased poleward movement in the eastern 
BPeaeitfic. This in turn appears to weaken the influence of 
Meeoemidad Pacific trough increasing the anticyclonic circu- 
lation over the mid Pacific and strengthening the easter- 
[eee into a continuous band across the Pacific. This 
sequence of events parallels Bjerknes (1969) hypothesis. 

We have shown correlations do exist among parameters within 
the eastern Pacific and the western Pacific regions and 
also between the two regions. Since the best correlations 
occurred between parameters calculated over the same geo- 
graphical area we feel more attention given to area 
selection would produce more meaningful results than by 


just sampling each parameter in each region. 
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Figure 19. Sea surface temperature anomalies 1974-75. 
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